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ABSTRACT 



An investigation of conventional and broken-symmetry scaling, 

and their effect on the ratio of the neutron-proton electric form 

factors, /G^ , was carried out using measurements from elastic 
’ En Ep 

electr on-deut er on and electron-proton scattering in the range of 

2 2 -2 

momentum transfers, q , of 0 . 10 ^q^ 0.80 F , and electron scat- 
tering angles of 4 5° to 120°. Within experimental errors of 
available data in this range of momentum transfers and angles, it 
was concluded that both forms of scaling as applied to the deuter on 

form factor have negigible effect on the ratio G^ /G^_ , and the 

En Ep 

application of either scaling law to the proton form factor results 

in nearly identical corrections to this ratio, thus allowing no 

distinction between them, and the agreement of the experimental 

2 

data with the neutron-electron interaction slope at q ^ 0 is not 
influenced by the choice of scaling lav/. A conclusive test for 
the scaling laws would require an improvement in experimental 
accuracy by at least an order of magnitude, a technically nearly 
impossible requirement. 
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I. 



INTRODUCTION 



The elastic scattering of electrons by protons is described in 
the first Born approximation by the Rosenbluth equation [lj. In 
this approximation a single virtual photon is exchanged between the 
electron and the proton and the structure of the proton is repre- 
sented by two electromagnetic form factors which are functions of 
2 

q , the square of the four -momentum transferred by the virtual 

photon [2]. Janssens, et.al., L3] made absolute measurements of 

the elastic electron-proton cross sections with an accuracy of 

about 4 % in the range from 4.0 to 30.0 F ^(1 F = 1 Fermi = 10 ^ra) . 

These measurements verified the Rosenbluth equation for values of 
2 -2 

q up to 22.0 F , and determined the form factors of the proton 

to a precision exceeding that obtained in previous experiments. It 

was also found that the charge and magnetic form factors of the 

2 

proton have the same dependence on q within the accuracy with 
which they were determined. 

Hundreds of other elastic electron-nucleon scattering measure- 

2 

ments are also available, extending over a very wide range of q 
and angle 14 ] . Empirically, the data seem to exhibit remarkable 
simplicity, the electric and magnetic form factors of nucleons 
obeying the "scaling law"; 



G (q 2 ) = 

Ep V 4 ' p 






= G 



SL’ 



Sn^) = 



0 



( 1 ) 



with the 



2 

q 



dependence being well described by the "dipole fit": 
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1 



( 2 ) 



J SL Dipole 



(1 



2 2 
q /m^) 



2 -2 

where = 90.25 F . It has long been clear that the electron- 
proton data show systematic oscillations about the dipole fit 
[ 556 ], but the amplitude seems to vary with laboratory and ( 2 ) is 
regarded as a remarkably good phenomenological fit. 

Following the same general procedure as Rosenbluth did with 
the proton, Jankus [ 7 ] derived the cross section for the elastic 
electr on-deuter on (e-d) scattering from a point deuteron in the 
first Born approximation. His results differ qualitatively from 
the proton cross section in one important way, i.e., there are 
three form factors characterizing the deuteron while two suffice 
for the proton. This is in agreement with the work of Glaser and 
Jaksic [8]. Their study showed that the cross section for the 
scattering of a relativistic electron from a potential with spin J 
contains 2J + 1 form factors. The three deuteron form factors 
correspond to the charge, quadrupole, and magnetic moments in the 
static limit q 0 . Dr icker and Hand [ 9 ] experimentally verified 
for the deuteron an emper ical relation between the form factors 
which is analogous to that for the proton. The scaling law for the 
deuteron may be written; 



ED 



MD 

[i 1 
D 



M 

p.l = 1.71 = 7“ M* 



D 



M 



D 



(3) 



Based on these scaling laws, Stewart [lO] in a study of the 
neutron charge from factors at very low momentum transfers, found 
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that : 



(a) The Feshbach-Lomon wave functions [ill, together with 

relativistic corrections, removed an apparant discrepancy between 

2 

the neutron-electron interaction slope at q ^ 0 , and the slope 

given from values of (obtained by electron scattering) in 

2 -2 

the range 0.10 ^ q ^ 0 . 80 F 

(b) Within the relatively large errors propagated into G , the 

Bn 

Partovi wave functions [ 12 ] with relativistic corrections applied, 
are in disagreement with the neutr on- electr on interaction slope. 
The important conclusion from Stewart’s experimental results 
is that the neutron has a non-zero charge form factor, implying 
a charge distribution within the neutron. 

Schumacher [4] however, points out that the theoretical sit- 
uation is not satisfactory. The scaling law (1) led to the in- 
vention [l3] of the phenomenlogical symmetry now called SU(6) , 

w 

but a deeper interpretation appears to be lacking. Likewise, the 

dipole formula (2) corresponds to a double pole of mass 843 MeV 

instead of a combination of single poles at masses of realistic 

vector mesons. Although it was hoped that the oscillations around 

(2) might be explainable in dynamical terms, Goitein, Dunning and 

Wilson [ 5 ] showed that existing theories were not adequate to do 

so and that form factor dominance by vector mesons led to trouble 
2 

at high q . There has been little change in the situation they 
described, e.g., the best subsequent models still not being for- 
mulated such that an isospin reflection connects realistic de- 
scriptions of the proton and neutron [4]. Because of this awk- 
ward situation regarding suitable dynamical models for the form 
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factors, it is desirable to ask the data to provide as many 
answers as possible to questions which can be formulated without 
a committment to the restrictive details of particular models. 

Schumacher examines the symmetry content of (1) and its 
dynamical implications, and finds that the form of the deviations 
from (1) suggests a particular type of symmetry breaking, well 
known in the physics of strong interactions. From these consider- 
ations a broken symmetry theory of the form factors is formulated. 
Using the world’s e-p data, Schumacher obtains several measures 
of symmetry breaking, which imply an oscillating G whose sign 
and magnitude agree with the Krohn-Ringo slope [l4l. He further 

derives a set of scaling laws which relates the form factors G. 

Mp 

and G w to G_ and G^ . These relations are 
Mn Ep En 







G Mp 


— II. Cl 

"P^Ep 


+ Li G 
n En 


( 4 ) 






G Mn 


II 

■P 

a 

"d 


+ U G 
p En 




EP 


= g es 


+ g ev 


is the 


’’symmetr y-conser ving” part 


and 


ES 


> 

W 

O 

11 


is the 


’’symmetry- violat ing M part of the 


form 



factors. These conclusions are primarily model- independent and 
do not require pole dominance of the spectral functions. 

It is then, the purpose of this thesis to: 

(a) Test the two forms of the scaling laws using the data of 

r n 2 -2 

Bumiller and others L19J in the range 0.10 ^ q ^ 0.80 F , and 

(b) To introduce a method of finding the ratio G /G without 

En Ep 

using the scaling laws, thus finding information about Gp n . 
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II. THEORY 



A . GENERAL 

Jankus [ 7 ], found that for the deuteron 



dcr 

da 



= (^~ N ) (1 + 

\ ^ ' Mott \ 



2E . 2 © 
2 S 1P 2 



,-ir 



m d c 



I T7( 1 + 2 ( 1+J7 )tan 2 §) D^(q 2 ) 



2 , 2 8 22 , 2 

D c (q ) + 9 D Q (q ) 



(5) 



wher e 



D c ( q2 ) =J [ u2 ( r ) + w 2 (r)1 j o (T)dr 



( 6 ) 



is the charge structure factor , 



2. _ 6 n |~2~ 2 



D (q ) = 2 

q 



M 



2 c 

D J 



u ( r ) w ( r )• 



2 / \ 
w . ( . r ) 

W 



5 2 (T ) dr (7) 



is the quadrupole moment structure factor, and 



D M (q2) = i [u 2 (r)-% w 2 (r)] j o (T)dr + [u(r )w(r 



f8- 



j 2 (r)dx 



(8) 



accounts for the contribution of the intrinsic magnetic moments 
of the proton and the neutron to the scattering process, and 



^(q 2 )=|( w 2 (r> 

o 



j 0 (T) + j 2 (T) 



dr 



(9) 



is the magnetic contribution to the scattering process arising 
from the convection of charge in the deutron. Together (8) and 
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(9) are 



V q } = h 



M F 

(U. + p. )2 D m + D 
p n M M 



The static limits (q 0) are 



D c (q 2 ) 



D Q (q 2 ) 



M . 2. 
D M (q > 



2 

M d Q 



1 - — p 

2 D 



( 10 ) 



E . 2. 

°M (q > 



3 p 

2 D * 



The j and j^ are spherical Bessel functions, T = |q |r/2, Q is 
the quadrupole moment, the percentage D state, and the u(r) and 
w(r) are suitably chosen wave functions of the deuteron. 

Adler [l5l, shows that in the impulse approximation, the 
charge, quadrupole, and magnetic form factors are given by 
(along with their static limits): 

* 1 ( 11 ) 



G_ s (G_ + G_ ) D — 

C v Ep En C 2_ 



0 



and 



= ( G_ + G_ ) D^. 
Q v Ep En ' Q 



■* m d q 



( 12 ) 



and 

f m ^ n 

G = (G + G ) D + (G + G )2D > — — U (13) 

M v Ep En ' M V >lp Mn ' M M D ' ° ' 

P 

where the D f s are given in equations (6) through (9)- The step 
from a point deuteron to a finite size deuteron is thus carried 
out by the introduction of the free nucleon form factors. This 
step is justified by the impulse approximation. 
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Rewrite equation (5) as 



(£>&) 



V”"/ Mott 



= A(q 2 ) + B(q 2 ) ta n 2 § 



(5a) 



where 



. 2 . 2 8 2 2 2 2 
A(q ) = G C + 9 V % + 3 V G M 



(5b) 



= a. 



, 2 2 
Z,M d c 



and 



B(q 2 ) 



4 2 

3* G m 



(5c) 



where the approximation 1 + r) 1 has been used. Equation (13) 
can be simplified considerably by scaling, i.e., 



and 



G = p T G , p, 1 = 1 .71 
MD D CD’ D / 



G G 

Mp _ Mn _ 

Pp [in ^Ep 



For the case of the deuteron 



G MD ^ G Ep + G En^ M, D D c* 



( i4) 



Then 



2 , 2r 2,, 2 ,2 4 ,2 2 9 8 2 2 -i 

G 0 = <<V G En> CD C< 1+ 3’’^ * tan 5> + ? n °Q ’ 



or 



2 2 2 

G = (G + G ) F* 

D v Ep En ' D 

,2 



(15) 



where F^ is called the deuteron structure factor and is the same 
as the result of Jankus for point nucleons in the deuteron. From 

(15) 



^ = F . 4 ♦ is' 

G_ D \ G 

Ep \ Ep 



(16) 



Again consider equation (5). If the quadrupole term is 
neglected (a contribution of order 10 ^ or less) and 1 + Tj ^ 1 
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is used, (5) becomes 




Mott 



g d = g e + f 77 (1 + 2 ta " 2 |) g m 



, 2 2 2 , 4 2 2 G 

< G E + 3 T,G M tan 2 



( 17 ) 



From this equation, using 



(£)'(&) 



Mott 



as the ordinate 



2 

and tan 9/2 as the abcissa, a linear plot can be made ( "Rosenblut h 
2 2 

plot”). Then is given by the slope, and can be found from 

the intercept. This method, however, presumes that absolute mea- 
surements of the deuter on cross section are available. Since 
Stewarts measurements are relative, available data must be an- 
alyzed in an attempt to make the experimental cross sections 
absolute in a first order sense. 

This task can be accomplished by considering the work of 
de Vries [l6] and Buchanan C 4]. De Vries collected all available 

proton data, and made a fit to these data using adjustable para- 

2 

meters so that a minimum was obtained in a X fit. Buchanan 
collected the world's proton, neutron, and deuter on data, and 
updated them. He checked the de Vries b T fit using the reanalyzed 
data, and concluded that this fit is the best to describe the 
proton form factors. 

Stewart [lO] measured 



R = 



da 

d-TL 



exp 



/ 



D 



dCT 

d.n 



exp 

P 



(18) 



If it is assumed that the de Vries b' fit gives G absolutely, 

P 



then 



SfT- .* 



abs 



da 

d -a 



i Mott 
'P 



(19) 
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If the ratio 



N 





exp 

P 



( 20 ) 



is formed, N can be considered to be the factor required to nor- 
malize the experimental data to absolute data in a first order 
sense. Then, since the ratio R was measured 




theor 

D 




( 21 ) 



can be considered to be absolute in the first order. The sub- 
script theoretical is used to distinguish absolute in the first 
order from absolute. 

2 

Thus Rosenbluth plots for the deuteron may be made, and G^ 

2 

and G can be extracted. If the scaling law holds, then within 
MD 

experimental error 



MD 



|JL 1 G 
D ED 



1 . 



B. SCALING APPLIED TO MAGNETIC CORRECTIONS 

2 2 

To eliminate the magnetic contribution to G^ and G p , the 
magnetic contribution is given by equation (15) as ■ 



C D = fVD ^ 1 + 2 ta " 2 f> 

so that 2 

2 _ G P eX P 

G ed (i+c D ) 

with 



2 /da 

G D e*P = ( — 



exp 



d si 



/ 



da 

d _n_ 



Mott 



D \ / D 

where the quadrupole terms has been neglected. 



( 22 ) 



( 23 ) 
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For the proton, the magnetic contribution is found by con- 



sider ing 



= _§£_ * (1+2 tan 2 f) . 



H77 



P 



Scaling assumes G.._ = M- G,^_ , where p = 2 . 973 > so, 

Mp p Ep p 



g p = 1^- i 1 + Vp (1 + 2 ta " 2 f»- 

p 



The term 



2 2 0 
C = 7? P (1 + 2 tan -) 
P P P 2 



is then the magnetic correction term for the proton. 
_ G 2 (l + 77 ) 

G 2 = _e_ lei 

Ep 1 + C 



wher e 



G 2 = 

P 



d_0_ 

d /l 



P 

exp 



/ 



da 

d si 



, Mott 



From his experiment, Stewart Cio] calculated 



( 24 ) 



(24a) 



( 25 ) 

Thus 

(26) 




1 g ed 
f d g ep 



1 



and evaluates from his data 




E K 



D 



F K 



( 27 ) 



where the E and F contain radiative corrections, Mott cross 
sections, target parameters and experimental counting rates. 

Of interest here is the ratio K^/K^, which represents the 
magnetic correction terras, it is given by 



K n 1 + C 

- E 

K ( 1+C ) ( 1+T) ) 

p v D' v ' p ' 
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where and were defined by (22) and (25). This ratio contains 
the scaling law assumptions. 

Next consider the scaling law proposed by Schumacher C4J • 
Specifically 



G “ JJL G + LL G 
Mp p Ep n En 



(29) 



G = M> G + M- G 
Mn n Ep p En . 

For the deuteron, consider similiar scaling laws 



g md 


It 

■P 

a- 


G + 

ED 


% 


G En 


(30) 


G Mn 


= M* 

n 


G + 
ED 


\L f 
D 


G En . 





Then equation (17) becomes, neglecting terms containing G , 

En 



Ep 



G„ +G, 



2 ['. 


2 


f 2 


(i 


„ 2 


© 


c L ] 


L + — 71 

3 ; D 


ji. 1 
D 


+ 2 tan 


2 


f 4 


T? M* M* T 
D n D 


■ (i 




0 2 © 






D 


+ 


2 tan — 


)) 


V'' 


C 








/ J 



(31) 



(32) 



Now define 



s = /l v n ■ + i V Gl5n v.A(i + 
o (.sVd d c »d)[ 



2 tan 2 | 



(32) 



Similar ily, for the proton equation (2^) becomes, with the use 



of equation (29), 

1 + (n U 2 + 2t) u u, )/l + 2 tan 2 % 

p p p P nG Ep A 2 



.2 _ 1E£_ 

p 1+ n, 



(33) 



and define 



c s = 
p 



(Vp * 



217 M* M- _ 

p p n G 



En 

Ep 



1 + 2 ta 



2 9 \ 
n 2 ) * 



(34) 



Thus when using Schumacher’s scaling, the magnetic correction 



ratio becomes 



+ C A 



(i+cp(i-n? p ) 



(35) 
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It should be pointed out that equations (32), ( 3^ ) , and thus (35) 

contain the form factors and , both of which are sought 

En Ep * 

from experiment. To use the readily available values from the 

de Vries b’ fit would be combining the dipole fit with the broken 

symmetry theories, however these are needed for comparison of the 

scaling laws. If the de Vries b 1 fit values of G^ and G^ are 

Ep En 

2 

incorrect at low q , their influence on equation (35) is quite 
small, so they will be used here. 

A comparison of equations (28) and (35) gives a measure of the 
effect of the two scaling laws on the data of Stewart. These re- 
sults are shown in Table VI. 



C. DERIVATION OF AN EQUATION FOR G^ /G^ WITHOUT SCALING 

En Ep 

This section has a threefold objective: 

(a) Derivation of an expression for G^ /G r _ without the use 

' En Ep 

of scaling laws, 

(b) Testing and comparison of the values of G^/G obtained 

2 -2 

from the equation derived in (a) above, for 0.1 ^ q ^ 0.8 F , 
using de Vries b 1 fit, or proton scaling wherever required. 

(c) Testing the effect of Schumacher’s scaling on G^/G . 

To derive an expression for G^ /G^ without the use of 

En Ep 

scaling laws, use is made of equations previously introduced, 
specifically 



and 



G _ = ( G^ + G^ ) D_ 
C v Ep En 7 C 



F M 

G==(G + G ) D + (G +G)2D 

M K Ep En 7 M v Mp Mn ; M 



M E 

where D^, D^, and D are defined in equations (6), (8), 

respectively, and 

2 2 2 2 0 

G d = A( c[ ) + B( qf ) tan' f 



( 11 ) 

(13) 

and (9) 



17 



(5a) 



2 2 2 2 

a ( q ) = g cd + i% g md 



= ^D G MD 



By the above equations, 



2 2 2 

G = (G + G ) D 
D ' Ep En ' C 

1 2 

+ - 71 T( G + G )D E + (G + G . ) 2 D^ f(9) ( 36 ) 

3 | ' Ep En ' M v Mp Mn y M J K y 



whet e 



f(9) = 1+2 tan 2 |- . 



For the proton, use equation (2 6 ) 

_2 



G 2 = + 

P 1+ Hr, 



f ( e > G Mp 



and recall from equation (l8) 



R = 



da 

d/i 



/ 



^ ( HL 



/ 



D 

2 ,2 . 



\ d.n , 

\ / P 



exT' 



2 2 

By taking the ratio G^/G^ it is possible to rearrange terms to 



find 

<4 . 




2 2 2 r F m ^ 

_ ^ G Ep +G En^ D C + 3 7 ? L^ G Ep +G En^ D M +2 ^ G Mp +G Mn^ °M f(9) 



-^P- .. + ^IL_ f(Q\ Q 2 

!+»] i+r] Mp 



Again rearranging terms, end up with, 
Mott 



R 



da 

dA 




Mott 

D 



H®1 + 
l *» P liT > P 



.2 /G, 




= D 2 (^+ 1 | F( 9 ) 



( 37 ) 



| 7? J D E“ /I^ + 2D B 2 rW\ +4 d H d M 

3 Dj M \ G_ I M G„ j M M 
L V Ep ) \ Ep / 
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2 






G +G 
My> M n 

G Ep 




wher e 

F( Q) = f (0)" 1 . 

This equation is free of approximations and does not contain 
scaling . 

To put equation (37) in to a form free of magnetic contr i- 

2 

butionsj a resort to graphical methods at constant q was made. 
The equation can be written in the form 



Y = mX + b 



wher e 



Y = R 



'do\ Mott 

- d Vn 



P(9) 

1 +J7 



1-^ 



< 



m = D 



En 



Ep 



+ 1 



F(9) = (1 + 2 tan 2 %) 



-1 



,2 G 



b = T 7 ? 



3 'D 



d; 



M 



En 

"Ep 



,2 G 



+ 2D 



M G 



En 

Ep 



G . +G. 



+ Z,d E d M _En _ M P _ M n 

M m(g epA g e p ) 



2 / G . +G 



♦ d b ~* Ad e d h ( * 4d m 7 ^p_^SL 

M M M l G ] M ^ Gep 




For constant q , X and Y are the only variable terms, Y = f(R,9) 

2 2 

and X = f(9). The assumption was made that the error in G.^/G^ 
x ' Mp Ep 

2 -2 

for 0.1 £ q ^ 0.8 F has a negligible effect when compared to 

2 2 

the error in R. It should be noted. here that G. /G^ has to be 

Mp Ep 

obtained either from the scaling laws of the proton, or fron the 



de Vries b f fit; a later discussion shows that it makes little 
difference which. 
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Using a least squares fit to the data for the above equations, 
the graphical methods were found to be relatively inaccurate. 

The remedy to the above difficulty was to turn to an analytical 
method where the points were averaged using a weighted average 
method. Equation (37) can be rearranged to give 



/dO_\ MOtt 

[dn ) 



'd a\ Mott 



F( 9 ) + 



i+n. 



_ r 2 

P .. 

Ep 



= |’!sd\Y d 2 f(9) * 2„ / 
1 %/ 1 c 3 D M 



a aj( 2D 2 F(e) ,| % ^*fn D D^D«sy 



(38) 



wher e 



- f D C F ( 9 ) + I D M 2+ ! D M °M s + ! ^ 



G Mn +G M n 

s = 



Ep 



The assumption was made that every quantity can be taken to be 

exact (although the D ! s are based on theoretical models), except 

2 2 2 2 
R, G %4 / G_ , S, and consequently the unknown, /G^ 

Mp Ep n En Ep 

Up to this point equation (38) is exact and contains no 

2 2 

scaling, however, G w /G_ , and S are as yet undetermined. There 

Mp Ep 

2 2 

are two methods available to get G /G^ since absolute data are 

Mp Ep 

not available. The methods are from the de Vries b T fit and the 
scaling laws. 

If the de Vries b f fit is used, 




S =( -MR — M£ \ de Vries . 

% 
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If scaling is used, 



4 P 



= y* 



S = P + P 
P n 



Now equation (38) i s put into a more workable quadratic form 

aX 2 + bX + c = 0 

where X = /G^ . By inspection 

En Ep 



a = Dc F( e > * f » D ^ 



2 4 8 F M 

b = 2D c F(0) + |t ?D D^ + - % D^D M S 



2 „E 2 8 „E „M _ 8 „M 2 _2 



c = D c F ( e> t - „ D d m + - „ D d m d m s + 1 d m s 



dO 
■R dA 



Mott 



JL 



f d0 

dst 



, Mott 



D 



n 



l+?7 + 1+17 

p p 



( 



% 

G Ep, 



The results of these two methods are summarized in Table VII. 

By applying Schumacher 's scaling to equation (38), a comparison 

of G^ /G^ can be made to that of the two methods discussed above. 

22 G M +G M 

In equation (38), the terms involving G x , /G_ , S = — £ — , 

Mp Ep °Ep 

2 

and S are the only ones affected. Let X = G^ /G . Then by 

En Ep 

application of equation (4)> it is seen that 

2 



^' 2 



Ep, 



= (P + P X) 
p n 



(39) 



and 



s = (M p + M- n )(l + X) . 
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By rearranging terms after substituting equation (39) into 



equation ( 38)5 the resulting equation is 



a f X + b 1 X + c' = 0 



(AO) 



wher e 



a' = ^ F(9) - fD D D® 




n 



P .2. 8 , » . 8 _M , .2 



T^ D m D m (* J + 7*n D M (** +*0 (AOa) 



1+7? n 3 DMM v p n 7 3 D M v p n ' 

P 



b' = 2D C D h « 



8 EM 

— V D D (|JL + [i ) 
3 D M M p 



+ 2 



8 M , . 2 

[S^D °M <W 



m 

-R 1 



Mott 

P_ 



J!e_ 



M- P 



do_\ Mott 1+r ? p P n 

<Wn 



(40b) 



and 



C = c = F( 9 ) | b D D* (b p * % ) 



8 M . , 2 

+ 7t) D I 1 + P ) 

3 D M p n' 




M£J + Jk. u 2 



l+n l+n p 
P P 



(40c) 
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III. REEVALUATION OF EXPERIMENTAL DATA 



A. DATA 

The data analyzed here were collected by Bumiller and others C 19 ! 

at the Naval Postgraduate School linear accelerator in the range 
2 -2 

0.1 ^ q ^ 0.4 F , and at the Mark III accelerator at Stanford 

2 -2 

University in the range 0.2 ^ q ^ 0.8 F . The experimental set 
up was basically the same in both laboratories. Whenever possible, 
overlapping sets of data were taken to check for any systematic 
errors in the equipment. No systematic deviations between the two 
sets of measurements were found and the data analysis was per- 
formed independent of the source. Table I shows the data used in 
these analyses, with the error listed. 



B. ROSENBLUTH PLOTS 



Equation (17) formed the basis for the Rosenbluth plots. This 



equation was broken up as follows for the plots: 



D 



D 



(17a) 



wher e 



2 0 

X = tan — 



'\theot = /do 

dsi I R(G p ) abs \ da 



Mott 

P 



and R is given by equation (l8). Wherever appropriate, weighted 
averages were used, and the corresponding error calculated. A 
weighted least squares fit was used to determine the slope and 
intercept, with the error in each. This method is described in 
Appendix C. 
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2 

Since these plots were made for a constant q , it is readily 

apparent that from the slope m, 

2 /„ 2 \ = 2 JL 

g md ’ 4 r? D 

and from the slope m and intercept b, 

g ed ( q ) = b ' 

Table III gives the results of these plots for the deuteron, 
and a sample plot is shown in Figure 1 . 

C. MAGNETIC CORRECTIONS 

The magnetic correction gives the number by which the experi- 
2 2 

mental G or is to be multiplied to remove the magnetic contri- 

2 2 

butions and yield G^ or G^ These corrections were computed 
^ Ep Ed ^ 

using equations (22), (25), (27), (32), (34), and (35). The 
terra contained in equation ( 32 ) is the charge structure factor 
and is given by equation (6), with the Feshbach -Lomom wave 
functions Cll] being used to describe the deuteron. Values of 
are listed in Table III. The magnetic corrections using both 
scaling laws, a comparison of the resulting corrections, and the 
effect of magnetic corrections using Schumacher scaling on G^/G^ 
as calculated by Stewart are shown in Table VI. The error was 
calculated by taking K^/K^ as the standard value. 



D. CALCULATION OF G /G^ WITHOUT SCALING 

En Ep 

Equation (38) is the basis for these calculations. The , 

E ‘ M . 

D^, and are respectively, the charge structure factor, magnetic 
contribution to the scattering process arising from the convection 
of charge in the deuteron, and the contribution of the intrinsic 
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magnetic moments of the proton and the neutron to the scattering 

process. These factors are given by equations (5), (8), and (9). 

As with the magnetic corrections, the Fes hbach -Lora ora wave functions 

were used. The values of these wave functions at the appropriate 
2 

q are also listed in Table III. 

2 2 

G /G is contained within equation (38) and is an undetermined 
Mp Ep 

quantity. Since absolute data are not available, this ratio may be 
found from the de Vries b' fit, or from the proton scaling law. 

Both methods of calculation were used, and it was found that they 
differ by a small percentage. 

Results of these calculations are listed in Table VII, and a 
comparison of G^/G computer using data from de Vries b T fit, 
and using proton scaling is shown in Figure 5. 



E. CALCULATION OF G^ /G^ USING SCHUMACHER SCALING 

En Ep 

These calculations were made using equation (40), with the 

coefficients given by equations (40a-c). As in the previous cal- 

E M 

culations the D , D , and D were calculated using the Feshbach- 
c M M 

Lomom wave functions as listed in Table III. 
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IV, PRESENTATION OF RESULTS, DISCUSSION, AND CONCLUSIONS 



A. PRESENTATION OF RESULTS, AND DISCUSSION 

The experimental data are presented in Table I. The errors 

quoted are attributed to counting and target thickness errors, 

and are the only errors that enter into a ratio experiment. 

Table II shows the Mott cross sections for point proton and 

2 2 

deuteron in the laboratory system; G^(q") from the b 1 fit of 
de Vries is the absolute form factor for the proton; the theo- 
retical cross sections as described in Section IIA; and the nor- 
malization factor N required to make the experimental deuteron 

cross section absolute in the first order sense. 

E M 

Shown in Table III are the quantities D , D w , D w as described 

^ c 5 M M 

in Sections IIA, computed from the Feshbach-Lomom wave functions 

without relativistic corrections. 

Table IV lists the magnetic and electric form factors of the 

deuteron as calculated from the Rosenbluth plots. Figures 2 and 

3 illustrate these results. It is seen from Figure 2 that the 

electric form factors approach the static limit as predicted by 

the theory. From Figure 3 it is seen that the magnetic form 

factors are to large by a factor of three. This is due in part 

to the fact that G is a function of the slope of the linear plot 

and is very senstive to experimental errors. A weighted least 

squares plot was used in an attempt to minimize the error. Also 

G 

shown in Table IV is the ratio — ; — ~ — . If scaling is to hold, 

^ d g ed 

this ratio must be unity, within experimental error. It is seen 
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however, that within experimental error — 



MD 



u, T G 
D ED 



3, indicating 



that scaling does not hold. Since the values of G are about 

MD 

three times larger than predicted, it must be held that this test 
is inconclusive. 

To investigate the effect of magnetic contributions to the 

electric form factor, and get information as to the importance 

2 

of scaling within the range of q of interest, the Rosenbluth 
plots were considered. From these plots it was seen that 

< 4 = 

where b is the Y intercept, and m is the slope of the linear plot, 



Specifically , 



C ED = b - ^ m 



where m contains the magnetic contribution. Define 

*2 

G = b 
ED 

* 

where G _ is the electric form factor without magnetic contri- 

but ions. Table V shows the results of these calculations. The 

result is that a very small error is introduced if the magnetic 

contributions are disregarded, thus concluding that scaling has 

2 

little effect in the range of q considered here. 

Stewart LlO] was interested in the ratio G /G^_ , from which 

ED Ep 

he finds G^ / . It was shown in Section IIB, that the ratio 
En Ep 5 

G ED /G Ep contains the magnetic correction ratio, which in turn 

contains the scaling laws. Using both scaling laws, a comparison 

was made of the two magnetic correction ratios. Table IV contains 

these results, showing that the difference in the corrections is 

2 

small, but tend to increase with q . 

Table VII shows four sets of values of G^ /G^ as obtained: 

En Ep 

(a) Without using scaling, but using the assumption of the de Vries 
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b T fit; (b) By using proton scaling in the form G x = M> : 
v 1 Mp p Ep 

(c) By Using Schumacher f s scaling law for the proton and; (d) As 

calculated by Stewart Clo] . In each case the Feshbach-Lomom wave 

functions were used, and relativistic corrections described by 

Stewart were applied. Each of the columns shows the slopes 
2 

d(G^ n /G^.^)/dcj , calculated by using a weighted least squares fit. 
In each case the slopes are in good agreement. 



B. CONCLUSIONS 

Based upon the experimental data available in the range of 
momentum transfers of interest here, the folloiving conclusions 
can be made: 

(a) Both forms of scaling as applied to the deuteron form 
factor have negigible effect on the ratio G^/G^. 

(b) The application of either scaling law to the proton form 
factor results in nearly identical corrections to the ratio 

G^ /G^ , thus allowing no distinction between them. 

En Ep 

(c) The agreement of the experimental data with the neutron- 

2 

electron interaction slope at q *** 0 is not influenced by the 
choice of the scaling law. 

Although both forms of scaling give satisfactory results in 
the momentum transfer range investigated here, the Schumacher 
scaling appears to be on a better theoretical basis, and should 
be tested at higher momentum transfers. A conclusive test for 
the scaling laws would require an improvement in experimental 
accuracy by at least an order of magnitude, a technically nearly 
impossible requirement. 
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The interpretation of measurements at higher momentum transfers., 
where the difference between the scaling laws would be more pro- 
nounced is hampered by increased difficulties in the theoretical 
prediction of deuteron wave functions. 

In agreement with Stewart, an important qualitative conclusion 
from these analyses is that the neutron has a non-zero charge form 
factor, not affected by the choice of scaling laws, implying a 
charge distribution within the neutron. The neutron will then have 
a charge structure with the outermost part of the charge distri- 
bution being negative in sign. 
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APPENDIX A. TABLES 



TABLE I. EXPERIMENTAL DATA 
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MOTT CROSS SECTIONS, G 2 ( h' ) , f--^) the ° r , N, ( theor 
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Weighted averages are used wherever applicable 



TABLE III 



M E 

D , D w AND D„ COMPUTED FROM 
C M M 

FESHBACH-LOMOM WAVE FUNCTIONS 



2 

q 

F - 2 


D c 


C M 


"5 


0.1 


0.9397 


0.8702 


0.0729 


0.2 


0.8866 


0.8214 


0.0717 


0.3 


0.8396 


0.7783 


0.0705 


0.4 


0.7974 


0.7398 


0.0693 


0.5 


0.7592 


0.7051 


0.0682 


0.6 


0.7245 


0.6736 


0.0672 


0.8 


0.6636 


0.6186 


0.0651 
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DEUTERON MAGNETIC AND ELECTRIC FORM FACTORS FROM ROSENBLUTH PLOTS 
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TABLE V 



EFFECT OF G_ >r . ON G„^ BASED ON ROSENBLUTH PLOT 
MD ED 



2 

q 

F - 2 


g ed 


2 

3 ” G md 


* 

g ed 


Error 

% 


0.1 


0.8919 


0.0294 


0.9083 


1.84 


0.2 


0.8536 


0.0110 


0.8601 


0.761 


0.3 


0.7986 


0.0074 


0.8032 


0.576 


0.4 


0.7775 


-0.0065 


0.7733 


o. 54 o 
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The magnetic correction gives the number by which the G or G is to be multiplied to 
remove the magnetic contributions and yield C‘l or G^r^ 

HD iCD 



COMPUTED WITHOUT SCALING, WITH PROTON SCALING, WITH SCHUMACHER SCALING, AND AS COMPUTED 
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Figure 
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Figure 
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Figure 
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Figure 
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APPENDIX C. STATISTICAL TREATMENT OF DATA 



A. WEIGHTED METHOD OF LEAST SQUARES 

The principle of least squares can be stated as follows: The 

most probable value of a quantity is obtained from a set of mea- 
surements by choosing the value which will minimize the sum of 
the squares of the deviations of these measurements. Specifically, 
if two variables X and Y are known to be related by a linear 
equation of the form Y = mX + b, where m is the slope of the line 
and b its Y intercept, and if a series of N observations ( x ^>y^) 
are made in which random errors occur only in the y^ measurements, 
and the observations are not all samples of the same parent dis- 
tribution, then the deviations are weighted inversly as the 
variances of their parent distributions, and Young L17J shows 
that the quexntity to minimize is 



£ w.d 2 = Z -i 
1 1 2 



= E w . (mx . +b-y . ) 
1 ' i i' 



2 



(4i) 



i 



wher e 



1 



w. oc 




1 




is the variance in y^ Ll8_). Upon minimizing equation (4l) 



with respect to b and m, the normal equations are 



mE w.x. + b E w. = E w.y. 

li i li 



( 42 ) 



mE-w.x. + bE w.x. = E w.x.y. 

li ii ill 



Solving for m and b, it is found that 



m = M/D 
b = B/D 



( 43 ) 
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wher e 



M = 



E w . y . E w . 



1 1 



E w .x ,y . E w.x. 
111 11 



B = 



E w.x. 
i i 

v 2 

E w.x. 

i i 



E w . y . 

l l 



E w . x . y . 

ill 



and 



D = 



E w.x . 
i i 

v 2 
E w.x. 

l l 



E w . 

l 



E w.x. 
i i 



The variances in m and b can be calculated by use of the 



equation 



s n = £ f TT~ \ s " 

Q V ^y ± ) y;, 



( 44 ) 



2 2 
where is the variance of the mean of Q, S is the variance 
Q 

of the mean of y^, and so forth. Application of equation ( 44 ) to 

the present situation gives 

v 2 
E w.x. 

2 2 xi 

S, = S — - — 

b y D 



( 45 ) 



and 



wher e 



o o E w. 

2 _ 2 i 

S - S jr— 

m y D 



2 2 
S = w. S 

y i y. 



B. WEIGHTED AVERAGES 

If a quantity x is measured by N methods, obtaining values 
x^...x^, which are known to have standard deviations S^...S^ T , 
then a simple arithmetic average would not give the best value 
because it would make the x_^ equal in weight, while their errors 
are assumed to be different. A weighted average is then used as 



a means of determining the best value. The question is how to 
assign the correct value of the weight factor w^ to to give 
the minimum standard deviation in the average value x. Making 
a generalization of Beers [l8] it is seen that 

S w . x . 

— 11 

x = — ~ 

E w . 

1 




and 



wher e 




E 



1 



1 





C. PROPAGATION OF ERROR 

Young [ ] 7*1 shows that if a quantity Q is to be calculated 
from several observed quantities x^, and 

Q = f(x.) 

then 




where is the variance of the mean of Q, and S is the 
Q x i . 

variance of the mean of x^. Equation (47) can be shown to be 
true even if different numbers of observation are made on the x^. 

Application of equation (47) gives the variance of the mean 
for the appropriate functional form of Q. 
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